Vogel-Fulcher freezing in relaxor ferroelectrics 
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A physical mechanism for the freezing of polar nanoregions (PNRs) in relaxor ferroelectrics is 
presented. Assuming that the activation energy for the reorientation of a cluster of PNRs scales 
with the mean volume of the cluster, the characteristic relaxation time r is found to diverge as 
the cluster volume reaches the percolation limit. Applying the mean field theory of continuum 
percolation, the familiar Vogel-Fulcher equation for the temperature dependence of r is derived. 
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Relaxor ferroelectrics (relaxors) have long been at- 
tracting considerable attention in view of their unique 
physical properties^ A key feature of relaxors is the ap- 
pearance of a broad temperature and frequency depen- 
dent maximum of the dielectric permittivity e(T, to) as 
well as the absence of long range ferroelectric order in 
zero field at any temperature^ The high value of the 
quasistatic dielectric constant e in a broad temperature 
range around the peak temperature T m , and especially 
the giant piezoelectric effect^ make relaxors attractive 
for various technological applications. Another charac- 
teristic relaxor property is the extremely slow relaxation 
below T m , which signals the onset of relaxor freezing. 5 
The characteristic relaxation time diverges at the freez- 
ing temperature To according to the well known Vogel- 
Fulcher (VF) relation&l 



■mexp[U/k(T-T )], 



(1) 



where tq represents the inverse attempt frequency, U the 
activation energy, and Tq the VF or freezing temperature. 
For obvious reasons, Eq. (JTJ) only makes sense for T > Tq. 

The above empirical VF law (TJ) has been experimen- 
tally observed in a variety of other systems such as su- 
percooled organic liquids, spin glasses, polymers etc. Al- 
though many theoretical ideas about the origin of the 
VF law have been proposed in the pastj&^ i 10 ' 11 : 12 ! 13 ! 14 a 
derivation of Eq. (JTJ) at the mesoscopic level is still lack- 
ing. 

Much of the experimental and theoretical research on 
relaxors has been focused on compositionally disordered 
perovskites such as PbMg 1 / 3 Nb 2 /303 (PMN) and related 
compounds^ Since only disordered ferroelectric systems 
show relaxor behavior, it is clear that disorder and/or 
frustration are key factors in producing a relaxor state. 
In particular, charge fluctuations are responsible for the 
formation of polar nanoregions (PNRs) below the so- 
called Burns temperature Td^ In PMN, for example, 
one has Td ~ 600 K, while the dielectric maximum oc- 
curs at T m ~ 260 K in the quasistatic limit. The PNRs 
can be regarded as a network of randomly interacting 
dipolar entities with the corresponding statistical distri- 
butions of their size and dipolar strength. The resulting 
collective state is reminiscent of a magnetic spin glass — or 



rather an electric dipolar glass — and it has been shown 
that its static properties can be described in terms of 
the spherical random bond-random field (SRBRF) model 
of relaxor ferroelectrics i 16, 17 Experiments performed by 
neutron scattering 18,19 and NMR techniques 20 indicate 
that the average size of PNRs increases with decreas- 
ing temperature and saturates below T m , suggesting the 
possibility of a percolation-type transition into a frozen 
relaxor statei^ 

Relaxor dynamics is characterized by a broad distri- 
bution of relaxation times g(logr), and it appears that 
the freezing process is associated with the divergence of 
the longest relaxation time in g(\ogT)&22- A typical em- 
pirical method used to analyze the dielectric permittiv- 
ity e(w,T) is to consider the dielectric maximum tem- 
perature T = T m as a function of frequency u). One 
may then define a relaxation time r = 1/w, which is 
found to satisfy the above Vogel-Fulcher (VF) relation 
(JTJ. For example, in PMN oriented along [001], the ex- 
perimental parameter values determined in this manner 
are: r = 10~ 12 s, U/k = 911 K, and T = 217 K. 23 
On the other hand, the longest relaxation time T max was 
found to obey the VF relation ([1]) with similar parame- 
ter values, i.e., t = 4.3 x 10 -11 s, U/k = 970 K, and 
T = 224 KM 

The purpose of the present work is to present a simple 
physical mechanism for the VF-type relaxation process 
in relaxor ferroelectrics. We adopt a qualitative physical 
picture of the relaxor state below Td based on a network 
of PNRs embedded in a highly polarizable medium.— One 
can imagine that the medium consists of a number of fluc- 
tuating reorientable dipoles and/or small size polar clus- 
ters. The first possibility corresponds to dipolar glasses 
and order-disorder-type relaxor ferroelectrics, whereas 
the second one applies to displacive-type relaxors; how- 
ever, intermediate cases may as well exist. Each PNR will 
polarize the medium within a space region bounded by 
the correlation radius r c . As the temperature is lowered, 
r c is expected to increase; this process will continue until 
a number of PNRs of similar size start to merge, thus 
forming a polarization cluster. Eventually, freezing will 
occur due to the growth of both the size of PNRs and 
the correlations between them^ 3 - 

Let us now consider a PNR with a core radius r and 
assume that the polarization cloud associated with it can 



2 



be described by a power-law radial dependence, 

P{r) -P (ro/r) 3 ; r > r , (2) 

and P(r) = Pq for r < r$. The local electric field at some 
distance r is proportional to the Lorentz field 24 



(3) 



with a local field correction factor ip = 0(1). The 
field £ couples to dipolar fluctuations of the surrounding 
medium and at a distance r induces an electric dipole 
moment 



p = a£, 



(4) 



which is proportional to the polarizability a. For order- 
disorder relaxors one has a ~ /i 2 /fcT, where [i is the 
strength of the fluctuating dipole moment. Similarly, 



in the displacive case we can write a 



where e* is an effective charge, M the reduced mass, 



and 



aokT the frequency of a renormalized quasi 



harmonic mode, which becomes unstable at zero temper- 
ature. The associated change of the electrostatic energy 
is given by 



SE 



1 



~-a£ 2 . 



(5) 



If \8E\ > kT, the thermal fluctuations will be too weak 
to destroy the correlations between the dipole and the 
PNR and a bound state will exist. The correlation radius 
r c then corresponds to the limiting distance for which 
\SE\ ~ kT. Combining Eqs. ©, and (J3J we obtain 



T 



(6) 



where TJ = ip/iPo/S^/^keo for order-disorder relaxors, 
and similarly = tpe*Po /3^/2Maokeo for the displacive 
case. Formally, we can require that r c — ► ro as T — > Td, 
implying that TJ is of the order ~ T<j. 

It follows that the correlation radius of each PNR 
in each case scales with temperature t = T/T£ as 
r c ~ i^ 1 / 3 . Similarly, the correlation volume v c = 
47rr 3 /3 scales as ~ and the PNR dipole moment 
p c = (47r/3) Jq C P(r)r 2 dr as p c - Po ~ | log t|, where 
Po = -Po4to/3 is the core dipole moment. 

As the temperature is lowered, the PNRs will grow 
in size and gradually start forming a connected polar- 
ization cluster. If n is the concentration of PNRs, they 
occupy a volume fraction r) = 47mr 3 /3 = AirnT^/ST. 
When r\ reaches a threshold value r) p , the PNRs will 
merge into an infinite cluster — a familiar concept from 
the theory of percolation^ The temperature at which 
the infinite cluster appears is therefore T p = AnnT^ /3r/ p . 
The percolation threshold for hard spheres on a lat- 
tice in three dimensions (d — 3) is rj p ~ 0.35, and 
i] p ~ 0.294 for randomly overlapping spheres^ whereas 



rjp = 1/3 for random hard spheres in the effective medium 
approximation^ Experiments on PMNi£ show that the 
volume fraction of PNRs saturates at ~ 0.3 below T ~ 
15 K. 

The theory of continuum percolatior*21 predicts that 
the mean cluster volume v increases as r) approaches rj p 
according to a power law v ~ (r) p — f?) -7 , or explicitly 



V = v (l - Tj/rjp)' 



(7) 



where vq is the critical amplitude. Assuming that at high 
temperatures T ~ Td the mean cluster volume reduces to 
the average core volume, we obtain an estimate for the 
amplitude, i.e., v$ ~ 47rrQ/3. 

In the mean field case, which is applicable to systems 
with infinite effective dimensionality, one has 7=1. This 
may be a reasonable value for the present system, since 
the PNRs behave as a fully connected random-site net- 
work in view of the long range character of dipolar inter- 
actions. Using the above expressions for r\ and rj p we can 
rewrite Eq. as 



(1 - Tp/T) 



(8) 



As the volume of the cluster grows, it becomes in- 
creasingly difficult for its total dipole moment to change 
direction, and eventually the reorientation will be sup- 
pressed completely as the mean cluster size reaches the 
percolation limit. The relaxation time r for the reorien- 
tation of the dielectric polarization is usually described 
by the Arrhenius law r = ro exp(U/kT), where the ac- 
tivation energy U is determined by the potential bar- 
rier for the relaxing particle to jump out of the potential 
well. For magnetic clusters in disordered magnetic mate- 
rials, Neel2& suggested that the activation energy could 
be written as the product of mean cluster volume v and 
an anisotropy factor Q an , 



t = T enp(vQan/kT). 



(9) 



As already noted by several author s^i^ an analogous 
relation should be applicable to PNRs in relaxors. In- 
serting the mean cluster size v from Eq. © into Eq. ([9]) 
we immediately obtain 



r = T e^p[v Q an /k(T ~T p )}. 



(10) 



This result has precisely the form of the VF equation (JTJ 
with U — v Q an and T — T p . 

The anisotropy constant Q an in Neel's formula © has 
the dimensionality of energy density. Its magnitude can 
be estimated by assuming that the core radius ro is of 
the order, say, ~ 2 nm. Using the value U/k — 970 K for 
PMN2 2 - we thus find Q an - 2.3 x 10~ 3 eV/(nm) 3 . 

The VF relation pop for the relaxation of dielectric po- 
larization in relaxor ferroelectrics has been derived here 
on the basis of a plausible power-law model for the po- 
larization distribution within the PNR. The principal 
mechanism responsible for relaxor freezing appears to be 
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the growth and percolation of PNR clusters^ culminat- 
ing in the formation of an infinite cluster. Alternatively, 
the system could be described in terms of random nor- 
mal modes, which are determined by the eigenstates and 
eigenvalues of the random interactions between PNRsj 3 -^ 
however, the relaxation time entering the equations of 
motion for these modes must again be associated with the 
growth of the PNRs discussed above and should, there- 
fore, obey the same VF relation. Thus, according to our 
present model, the two physical pictures, namely, that of 
cluster growth and of the freezing of local modes^^ are 
essentially two parts of the same general scenario. 

It should be noted that the above results are indepen- 
dent of any specific spin glass-type static model such as 
the SRBRF model. Of course, the SRBRF model re- 
mains applicable to true static phenomena, for example, 
the temperature dependence of the quasistatic dielectric 
response as observed in a field-cooled (poled) sample. 

The divergence of r at T p does not imply that the mo- 
tion of all PNRs is completely frozen for T < T p . Namely, 
the complex dielectric permittivity e(w,T) remains finite 
below T p , indicating that some degrees of freedom are 
still active at low temperatures. In fact, as already noted 
above, only ~ 30% of the PNRs are involved in the for- 
mation of the infinite cluster. Smaller PNRs in the re- 
maining space may continue to undergo a similar process 
of growth and percolation, suggesting that at any tem- 
perature below T p an analogous freezing mechanism may 
apply. One can thus introduce a probability distribution 
of VF temperatures w(T p ) in the interval < T < T p 
which could, in principle, be transformed into a nontrivial 
distribution of relaxation times. A simple linear shape of 
w(Tp)2i then reproduces the main physical features, i.e., 
the finite value of e(u>,T) at temperatures below T p and 
a frequency dispersion of its real and imaginary parts, 
although the agreement with the observed behavior of 
e(tu,T) is — admittedly — only qualitative. 



At present, it is not yet clear whether the above ap- 
proach is applicable in a straightforward manner to anal- 
ogous systems such as random ferromagnets and spin 
glasses, supercooled organic liquids, structural glasses, 
etc. In magnetic systems, the magnetic dipolar inter- 
actions are much weaker than the short range exchange 
interactions, which are responsible for the formation of 
superparamagnetic clusters. However, in spin glasses the 
RKKY interactions actually fall off as ~ r~ 3 up to some 
oscillatory prefactor, and magnetic dipolar interactions 
are believed to be relevant in some spin glass systems. 32 
Thus the present model may serve as a convenient start- 
ing point for these systems. Meanwhile, in relaxor ferro- 
electric polymers, the above physical picture seems to be 
applicable without serious limitations^ For structural 
glasses and supercooled organic liquids, the possibility 
of elastic deformations and their interactions playing a 
leading role in the growth of correlated clusters should 
be explored. 

In conclusion, we have derived the Vogel-Fulcher (VF) 
relation in relaxor ferroelectrics by introducing a meso- 
scopic mechanism for the growth of PNRs. The basic 
idea is that thermodynamic stability of the polarization 
density with a power-law distribution P(r) ~ r -3 en- 
tails a temperature dependence of the correlation radius 
r c ~ T -1 / 3 . Thus a cluster of PNRs is formed on low- 
ering the temperature, and its mean volume v increases 
until the percolation limit is reached at some temperature 
T p . Using the mean field theory of continuum percolation 
we have shown that v ~ (1 — T p /T) _1 , and since accord- 
ing to Neel the activation energy U for the reorientation 
of the cluster polarization is proportional to v, the VF 
relation follows immediately. 
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